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Abstract. We study nonsymmetric tridiagonal operators acting in the Hilbert space 
£'^ and describe the spectrum and the resolvent set of such operators in terms of a 
■■M I continued fraction related to the resolvent. In this way we establish a connection 

C^ ■ between Pade approximants and spectral properties of nonsymmetric tridiagonal op- 

erators. 
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oo . 1- Introduction 

O 

CN ' In this paper we consider operators in the Hilbert space £^, with the following 

^-v ■ representation in an orthonormal basis of this space: 

5: ^/?o 70 ...^^ 

t^ : (1) A 



ai (3i 7i 
a2 /32 72 

/\ • where ctfc, /3fc, 7fc are complex numbers, ctfc 7^ 0, 7^ 7^ 0, /c G N. The operator A is 

j^ ■ defined for any finite vector x = xogo + xigi + . . . + Xng-n in the orthonormal basis 

{gn}o^, and its domain of definition D{A) is dense in i"^. 

The symmetric case, where ak = '^k-i and (3k real, can under an appropriate 
chosen basis be reduced to an infinite Jacobi matrix enabling a deeper examination 
of this spectral theory. The main tools of investigation in this case are the classical 
moment problem (including the theory of selfadjoint extensions of unbounded sym- 
metric operators), the theory of general orthogonal polynomials and the spectral 
theorem for selfadjoint operators ([1], [2], [9], [11]). This connection between the 
spectral theory and analysis is fruitful for various points of view. For example, the 
scattering problem for a Jacobi matrix can be treated on the basis of strong (or 
Szego type) asymptotic results for orthogonal polynomials ([4], [6], [10]). On the 
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other hand, the perturbation theory gives new results for orthogonal polynomials 

([12]). 

In the general, nonsymmetric, case we can not use the spectral theorem and 
orthogonal polynomials, but in return we will obtain results of use for Pade ap- 
proximants and the theory of continued fractions. We denote as usual by <7{A) the 
spectrum of the operator A, Q{A) = C\ (J (A) is the resolvent set and the resolvent 
is -R(A) = {XI — A)~^ for A G O. By {(7n}o° we denote the orthonormal basis in £^. 
The function 



(2) 



(j)iX) = {R{X)go.go) 



is analytic on the resolvent set Q{A). If the operator A is bounded (this is the case 
when ctfc, (3kj 7fc are bounded), then the function (j){X) is analytic for |A| > ||A||. For 
convenience we have changed the basis of representation of the operator A. The 
new basis {e^}^ is defined by e^ = gn/dm where dn = 7o7i • • -Tn-i, do = 1. In 
this basis the operator (1) has the following form: 

/ Po 1 

ai7o (3i 1 

a27i (32 1 



(3) 



A 



\ 



V 



/ 



If we set ttr 



(4) 



Oinlr, 



A 



(3m then we have in the basis {cn}: 



•••/ 



/bo 


1 








ai 


bi 


1 








a2 


62 


1 



V 



The Chebyshev algorithm (Algorithm 7.2.1 in [8, page 248]) applied to (/){X) at 
infinity gives us the following continued fraction (J-fraction): 



(5) 



ail 



«2| 



\X — bo |A — 61 |A — 62 



The numerators Pn and denominators Qn of the nth convergents for this fraction 
satisfy the three-term recurrence relation: 



(6) 



anUn-l + bnUn + Un+l = Xl/r 



with the initial conditions 



Q-i 
P-i 



0, 
1, 



Qo 
Po 



n 



0. 



0,1,2, 



The main problem considered in this paper is to describe the spectrum and the 
resolvent set of operators of type (1) in terms of the continued fraction (5) and 
the monic polynomials P„, Qn and in this way to establish a connection between 
Pade approximants and spectral properties of nonsymmetric tridiagonal operators. 
In our main result, theorem 1, we state the criterion for the resolvent set Q{A) 
in terms of the growth of the polynomials Qn{z)- The proof of the necessary 
condition is essentially based on the results of Kershaw and Demko, Moss and 
Smith about decay rates of inverses of banded matrices ([5]). This theorem gives us 
analytic properties of the polynomials Qn{X) and the related Pade approximants 



Tj //n 



J-U„ „„J- O/' /l^ -\TiT„ 



^■i-^ — ^ J-U:, 



J-U„ f„11. 
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Theorem 2. If A is bounded and 

< Ittfcl, It/cI < Ci < oo 

for some Ci, then for any X G 0(A) there is a subsequence of Fade approximants 
Tin which converges to (j){X) with a geometric rate. 

In section 2 we formulate and prove theorem 1 and related results. Section 3 will 
be devoted to an application of these results to the behavior of Pade approximants. 

2. Analysis of the spectrum and the resolvent set of the operator 

Suppose A has a representation (4) in an orthogonal basis {cn}, where ||en|| = hn- 
We assume that the operator A is bounded. In this case A is defined on the whole 
space i'^ and A G Q{A) if and only if the following conditions are satisfied: 

(1) Ker{XI-A) = {0}, 

(2) Vn, 3zn : {XI - A)zn = e^, 

(3) ||(A/ — A)~"'^x|| < C||x||, for all finite vectors x = xqCq + xici + ■ ■ ■ + XnCn 
and some constant C. 

Proposition 1. Xq is an eigenvalue of A if and only if 

oo 

(7) 5^|g„(Ao)p/.^<+oo. 

n=0 

Proof. We can formally write Ax = Xqx and obtain for x = xqCq + xici + ■ ■ ■ the 
following system 

boxo + xi = Aoa;o 
aiXo + biXi + X2 = Aoa^i 

a2Xi + b2X2 + X3 = XqX2 

If xo = then Xn = 0,\/n E N. If xq 7^ 0, then Xn = QniXo)xo and in this case 
xeeiSj:'^\Qn{Xo)\^hl<oo. D 

Proposition 2. For given X, the equation 

{XI - A)z = en 

has a solution in £^ for all n > if and only if there exists a complex number 7 
such that: 

00 

(8) 5^|g.(A)7-Pn(A)|2/i2<+oo. 

n=0 

Proof. First we find the formal solution z^^' of the equation 

r \ T A\ 



4 

Let ^(0) 

(9) 



A.I. APTEKAREV, V. KALIAGUINE AND W. VAN ASSCHE 




in the basis {e^} then 



hizi^o 

&2^2,0 



+ ^1,0 
+ Z2,G 
+ ^3,0 



A2i,o 

A22,0 



So Zn^Q satisfies the same recurrence relations as (5n(A) and Pn(A) for n > 1, only 
the first relation is different. This implies that Zn,o = QnWj — PnW for n > 1 and 
for any 7 the first relation is also satisfied: 6o7 + 7(A — bo) — I + I = A7. Thus the 
sequence (5n(A)7 — Pn(A) is a solution of (9) for any 7. We need a solution in £^. So 
Co is in the image of {XI — A) iff J^o^ IQn(A)7 — Pn(A)|^/i^ < +00 for some 7. Note 
that if A is not an eigenvalue of A then the convergence of this series is possible only 
for one 7. For the solution of the equation {XI — A)z = ei we note that {XI — A)eQ = 
Xeo — Acq = (A — 6o)eo — aici and then — aici = {XI — A){eQ — (A — 6o)-2*-*^'*)- So 
z'^ "> = — (eo — (A — 6o)-2 )/fti is a solution of {XI — A)^ = ei. This is an element 
of £^ if z^^^ G £^. In the same way we find z^"^^ = — (ei — (A — bi)z^^^ + z^^^)/a2 for 
the solution of {XI — A)z = 62 and so on. D 

Theorem 1. Suppose A is bounded, X is not an eigenvalue of A and (8) is satisfied 
for some 7. Then X G 0(A) if and only if there are constants C > and < q < 1 
such that 



(10) 



(11) 



, Qn{X)rm{X) hn 
G.IG.2 • • • O-n hj 

,Qm{X)rn{X) hr. 



aia2 



■ Ojri, ilr 



<Cq' 



<Cq' 



n < m. 



n> m. 



where rra{X) := Qm{X)'^ - Pm{X). 

Proof. From the propositions 1 and 2 it follows that the inverse operator B = 
{XI — A)~^ can be defined on the basis vectors by Bcn = z^"^' . The matrix of B in 
the basis {cn} is of the form: 



(12) 



Indeed, 
column 



fro 
r2 

^3 

V : 



ri 
a-i 

Qi 

■^ ai 



r2 
aia2 



ra 



^ ^ aia 
Q2 

rs 



r2 



aia2 
Q2 



aia2 



a\a2az 

0,—!-3— 
^ J^ aia2a3 

O-?—^ 

^^ aia2a3 
13 



Qs- 



a2as 



\ 



I 



for the first column there is no problem, since it contains z^^'. The nth 
contains z^'^\ For z^'^^ = — (eo — (A — bQ)z'^^^) j a\ we have 



1 



(13) 



^0,1 = 
^1,1 = 

Z2A = 



(A - bo)ro _ n 
ai 



-(A-6o)^i,o 

ai 
-(A-5o)^2,o 

ai 



ai 

Qi 



r2- 



Ol 
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This is the second column of B. Suppose we have calculated z^'^^ as indicated and 
we have to find z^'^'^^' . First note that 

ftfc+i 
Then 

1 , ,. , X rfc rk-i ■, Tk+i 



ttk+i aia2---ak aia2- ■ ■ ak-i aia2---ak+i 
^i,fc+i = (-(A - bk) \ ) = Qi 



ttk+i aia2---ak aia2---ak-i aia2---ak+i 



^fc,fc+i = (1 - (A - bk) \ ) = Qh 



ttk+i aia2---ak aia2 ■ ■ ■ ak-i a,ia2---ak+i 

We used here the relations Qk-ifk — Qkfk-i = —aia2 ■ ■ • dk-i deduced from the 
recurrence relations for Qk and r^. Next, 

1 . ., , N Qk^k+i , Qk-iTk+1 N Qk+1 

Zk+i,k+i = (-(A - bk) \ J = Tk+i- 



a-k+i aia2---ak aia2---ak-i aia2---ak+i 

and so on for zi^k+i, I > k + 1. Thus the form of the matrix B is correct. Now we 
decompose B into two parts: B = Bi + B2, where Bi and B2 are lower and upper 
triangular matrices which the same diagonal elements {Bi)i^i = (-62)1,1 = ^Bi,i- 
For the estimation of norms of matrices we use the formula 

||i^|| = sup sup \{y,Kx)\. 
\\^\\<i \\y\\<i 

We have for two finite vectors x = xqCo + xici + ■ ■ ■ + XnCn and y = yoeo + yiei + 

■ ■ ■ "r yn€-n- 

1- 



(y.Bix) =-bQ^QyoXohQ+ 



{bifiUiXo + -bi^iyixi)h-^ + 



1_ .. 



+ (bufiUnXQ + bn,iynXl H h -bn,iynXn)K- 



Here, we denote by bij the elements of the matrix B. Then 

.. n n n—j 

\{y,Bix)\ < -^\bk,k\\^k\\yk\hl + ^^\bk+j,k\\xk\\yk+j\hl+j 

k=0 3 = 1 k=0 

^ n n n-j ^ 

= - ^ \bk,k\\xk\\yk\hl + XlXl \bk+j,k^--^\\xkhk\\yk+jhk+j\ 



k=0 3 = 1 k=0 

<C{-+q + q^ + --- + q'')\\x\\ \\y\\ 

l-q-+^ 
< <-^ — :: \\x\\ \\y\\. 
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A similar estimation is applicable for the matrix i?2. The first part of the theorem 
is thus proved. 

To prove that for A G fl{A) the estimation (10) holds, one can use the known 
theorem on the decay rates of the inverse of banded matrices (see [5]) and the 
matrix representation of B = {XI — A)~^ obtained in the first part of our proof. D 

3. Application to Pade approximants 

In this section we consider some applications of theorem 1. First of all we note 
that theorem 1 together with propositions 1 and 2 gives a characterization of the 
resolvent set of A in terms of the polynomials Pn and Qn- On the other hand this 
theorem also gives a result on the convergence of the continued fractions (5). We 
recall that the fraction Pn/Qn is a diagonal type Pade approximant for (/)(X) at 
infinity (see, e.g.. Theorem 7.15 (B) in [8, page 250]). 

Corollary 2. If X E ^{A), then the remainder of Pade approximation in linear 
form tends to zero, i.e., 

lim[Q^(A)(/>(A)-P„(A)]/i, = 0. 

Proof. We have 2*^°-' = i?(A)eo and ^(A) = (i?(A)eo,eo) = zq^, with Hq = 1. On 
the other hand ^0,0 = iQo ~ -Po = 7? hence (/)(A) = 7 and the corollary follows from 
proposition 2. D 

Remark 1: In the case of a symmetric operator A with a representation (1) in 
some orthogonal basis, the polynomials 

Qn = Qn/dn, d^ = 7o7l " " " 7n-l 

are orthonormal and if A G 0(A) then qn{X)(j){X) — Pn{X) ^ as n ^ 00. In this 
case we have hn = l/\dn\. 

Corollary 3. // there exists a positive constant Ci such that \ak\j |7fc| < Ci then 
for all X G 0(A) 

limsup|Qn(A)/inr/">l. 

n— >oo 

Proof. Both Qn and r^ satisfy the same recurrence relation (6). Hence 

Qn-lfn — QnTn-l = —0102 ■ ■ ■ O^-l, 

or equivalent ly 

^ r ^n i^n ^ 7 ^n — 1 i^n—1 -. 

^n — 1'T'n—lT '^'^T ^n'^nT 7 J-- 

iT'TiO'l ' ' ' dn 'T'n-1 'I'n-lO'l ' ' ' O-n-l 'T'n 

From theorem 1 we get 

—^ <Cq', q<l, 
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SO the sequence Qnhn can not be majorized by a geometric sequence p" with p < 
1/q, i.e., for any p < 1/q and any positive constant C the inequahty 

\Qn{X)hn\ < Cp'' 

is not satisfied for an infinite number of indices n. For a subsequence A C N we 
thus have 

\QniX)hn\ > Cp'', n G A 

and consequently if we choose p such that 1/q > p > 1 we have 

hmsup IQuWhul^^"" > P > 1, 

n — >oo 

giving the required result. D 

Remark 2: In the symmetric case this corollary gives us the well known charac- 
terization of Q{A) in terms of the orthonormal polynomials q^: if A G fl{A) then 
limsup,_^|Q^(A)|i/->l([ll]). 

The combination of these two corollaries gives 

Theorem 2. if A is bounded and < |afc|, |7fc| < Ci < oo for some constant 
Ci, then for any X G 0(A) there exists a subsequence of the Pade approximants 
T^n = Pn{.X)/Qn{X) which convcrgcs with a geometric rate to (/){X). 

Proof. From corollary 1 we have 

lim Qnhn(l){X) - Pnhn = 0. 

n — >oo 

Let A C N be such that 

lim \Qn{X)hn\^/^ > 1, 

n — »oo 

then \Qn{X)hn\ > Cp"^ for some p > 1. Hence 



m ""-^'^ 



Qn{X) 



<c 

■ p 



which proves the theorem. D 



Remark 3: In the symmetric case this statement seems to be unnoticed. Indeed, 
suppose A is in the form (1) in an orthonormal basis, and consider the so-called 
asymptotically periodic case (or the corresponding limit-periodic continued frac- 
tion) where for some A^ we have 

lim ttfcAf+i = 5i, i = 0,1,2, ...A^- 1, 

k — »oo 

and similarly for (3kN+i and ■jkN+i^ then it is known that the poles of the Fade 
approximants Tin = Pn/Qn are essentially concentrated on a system of A^ closed 
intervals ([4], [7]). The convergence of n^ however depends on the behaviour of 
the so-called spurious poles. Nevertheless our theorem shows that for any A in the 
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Finally we note that from theorem 1 we get a connection between the fact that 
A G 0(A) and the convergence of Fade approximants 7r„(A) = Pn{\)/Qn{^)- The 
convergence of 7r^(A), without further requirements on the rate, does not imply that 
A G 0(A). A counterexample is given by the spectral measure d^{x) = Vl — x"^ dx 
on [-1, 1], for which P^(x) = 2-''+^U^_i{x) and Q„(x) = 2-'^t/^(x), where U^ix) 
are the Chebyshev polynomials of the second kind. At A = 1 G c^A) we have 
7r^(l) = 2n/{n + 1) which converges, but without geometric rate. The convergence 
of some subsequence with a geometric rate is also not sufficient to imply that 
A G 0(A). A counterexample is given by any symmetric measure on [—1, 1], because 
then at A = G a"(A) we always have P2n(A) = 0. 

The following question is thus of interest: does the convergence of the whole 
sequence 7r„(A) with a geometric rate imply that A G 0(A)? 
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